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Abstract. We present a new determination of the surface brightness of our Galaxy at the Solar Neighbourhood as observed from 
outside the Galaxy. We rely on various existing optical and infra-red surveys to obtain a multiwavelength estimate. On the one 
hand, scattered light does not contribute significantly to the surface brightness. On the other hand, optical and infrared integrated 
all-sky surveys (Pioneer 10/11 and COBE/DIRBE) show a systematically larger value than our synthetic local estimate based 
on Hipparcos data. This local estimate is also compatible with our Galactic simulations normalised at the Solar Neighbourhood 
and assuming an homogeneous stellar distribution. We interpret this disagreement as a signature of the presence of a local 
minimum of the stellar density compatible with Gould's belt. According to this result, the global luminosity of the Milky Way 
should follow the Tully-Fisher relation established for external galaxies. 
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1. Introduction 

Modelling of the multiwavelength continuum emission of high- 
z galaxies is essential to simulate realistically the formation and 
evolution of galaxies as will be observed with the next genera- 
tion of instruments. The data samples of galaxies are currently 
limited especially in the infrared and millimeter wavelengths, 
but will be significantly enlarged in the coming years. On the 
one hand, theoretical modelling based on physical principles 
actually relies on a large number of variables and requires well- 
sampled galaxy spectra to be constrained. On the other hand, 
empirical approaches (e.g. Guiderdoni et al. 1998; Dale et al. 
2001) minimise the number of parameters to reproduce the cur- 
rent observations. Both approaches need to be validated with 
well-sampled continuum spectra of galaxies. It is customary 
to use the "Solar Neighbourhood" as a zero-point for these 
models. Even though various multiwavelength observations are 
available for the Galaxy and the Solar Neighbourhood, most 
authors quote Mathis et al. (1981), who were investigating the 
Local Interstellar Radiation Field. Although this work is ro- 
bust to understand the absorption of photons by dust grains in 
the Local Neighbourhood, it did not intend to provide a spec- 
tral measurement of the surface brightness of our Galaxy at the 
Solar Neighbourhood as observed from outside the Galaxy. 
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In this paper, we compile various data sets to measure this 
quantity, which could be used as a zero-point for spectral mod- 
elling of external galaxies. We perform a cosecant averaging 
based on a plane parallel approximation for this measurement. 
In Sect. 2, we explain the basis for this approximation and its 
expected behaviour. We compare this approximation (integra- 
tion in a cone) with the method usually used for external galax- 
ies (integration in a cylinder) through simulations. In Sect. 3, 
we present a synthetic approach based on the Hipparcos data 
set and a stellar spectral flux library to estimate an optical 
spectrum of surface brightness of our Galaxy. We also present 
the two different types of integration (cone/cylinder), described 
here as corrected direct and indirect methods. In Sect. 4 and 
5, we present a compilation of various all-sky surveys with 
the cosecant averaging method. In the optical, we use the stel- 
lar catalogues Tycho-2 and USNO-A2 as well as the Pioneer 
10/11 integrated maps. In the infra-red, we project in cosecZ? 
the COBE/DIRBE survey. Last, we discuss all these estimates. 

2. Methods and simulations 

The purpose of this paper is to determine at different wave- 
lengths the surface brightness of our Galaxy at the Solar 
Neighbourhood as observed from outside the Galaxy. (Cf. 
Appendix B for a formal definition.) However, such an ob- 
servational measurement imposes strong constraints: (1) we 
want to compare different wavelengths (which could be sam- 
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pled differently); (2) we would prefer to avoid the Galactic disc 
area, which is affected by crowding (stellar catalogues); (3) we 
would like to sample a significant volume in order to reduce 
local fluctuations. 

The available all-sky surveys consist in stellar catalogues or 
integrated maps, at different wavelengths. In order to minimise 
the systematics due to this diversity, we adopt the same aver- 
aging procedure for each data set. The idea is to apply a cose- 
cant averaging on the data to estimate the surface brightness of 
the Galaxy. This method, also called integration in a "cone", 
satisfies the previous constraints, while approximating the in- 
tegration in a cylinder, usually used to compute surface bright- 
ness of external galaxies. In this section, we present the differ- 
ent aspects of this approximation with respect to the standard 
method. Firstly, we present the modelling performed for our 
Galaxy and the methods used to compute the surface bright- 
ness. We describe the "integration in a cylinder", which cor- 
responds to the standard method to compute a surface bright- 
ness. We explain the difficulties due to our position near the 
plane of the Galaxy, and the corrections required to estimate 
the surface brightness corresponding to a face-on observation. 
In parallel, we present the alternate method used throughout 
this paper, called here "integration in a cone". We discuss and 
compare the trends of both. Secondly, we simulate the cosecZ? 
projection that will be applied to stellar catalogues. 

2.1. Definitions 

We consider two extreme geometries for our Galaxy: (1) an 
exponential disc with scale length Iir = 2.5 kpc (Robin et al. 
1992) and (variable) scale height 1 hz(My), with a density nor- 
malised at the Solar neighbourhood; (2) a plane parallel model 
with a constant number density and a total thickness of 2h z = 
240 pc. In both cases, we assume that the Sun is located at a 
cylindrical Galactocentric radius r q = R B = 8 kpc, at z & = 15 pc 
above the Galactic plane (Humphreys & Larsen 1995) and that 
the total extension of the Galaxy with respect to the centre is 
Rgai = 14 kpc. The numbers n(r, b, I, My) of My stars per pc 3 
at the position (r, b, I) for the 2 geometries are as follows: 

n D.E X pon. (f . ^ ^ My) = n&(Mv) exp(ze/hz(Mv)) 

exp(R e /h R ) e\p(-r q /h R ) exp(-|z + z \/h z (My)) (1) 

with the cylindrical Galactocentric radius r q 
(r 2 + r 2 cos b 2 - 2r r cos b cos /) 
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the number density 
of My stars per pc 3 n (M v ) = <p(M v ) is the luminosity 
function 4>{M V ) at the Solar Neighbourhood, as defined in 
Appendix D; z = rsinb is the vertical position with respect to 
the Sun; 



n pp (r,b,l,Mv) = 



\n e (M v ) e\p(z /h z (M v )) if \z + z©| < h z , 







ifk + Zol > h z , 
(2) 



1 Similarly to Miller & Scalo (1979), we take h z = 80 pc for M v < 
0, h z = 250 pc for M v > 8, and a linear interpolation in between. We 
also consider a fixed scale height (hz = 120pc), compatible with the 
variable scale height modelling in Fig. 4. See also Appendix D. 
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Fig. 1. Simulations. Surface Brightness in V band of the Milky Way 
as a function of the cylinder radius R cy i for the models of Sect. 2. The 
cylinders corresponding to R cy i are centered on the Sun and perpen- 
dicular to the Galactic plane. The curves with symbols correspond to 
the surface brightness computed with a cylinder for the exponential 
disc or plane parallel geometry. We consider two different observers: 
one observing from the Solar Neighbourhood and one observing the 
Galaxy face-on. We compare an ideal galaxy without extinction with 
the real configuration with extinction. The arrows on the right-hand 
side correspond to the integration in a cone with the cosecant law ap- 
proximation, fitted for 1/ sin \b\ < 6. The caption for the line sizes (and 
types) of these horizontal lines are defined in Fig. 4. 



We account for extinction with a uniform dust layer following 
the prescription provided and discussed in Appendix E. 



2.1 .1 . Integration in a cylinder 

We define the average surface brightness vS v of our Galaxy at 
the Solar neighbourhood in the V band. 



v v Sy= [\ n(r,b,l,M v )f(A v (r,b))L(My 

Jr,b,l,M v 

r 2 cos bdrdbdl dMy] rq<Rcyl jnR 2 cyl 



(3) 



where L(My) is the luminosity 2 of a star of My magnitude, 
f(A v (r, b)) = Kr 4A "( r ' fo > is the extinction factor and R cyt is the 
radius of the cylinder used to compute the surface brightness 
(see Fig. 1). In the absence of extinction and for R cy i < 2kpc, 
the two geometries are equivalent and the surface brightness 
is nearly independent of the cylinder radius. With extinction, 
the surface brightness measured from the Solar Neighbourhood 
decreases with increasing cylinder radius. This is due to the ex- 
tinction affecting the \b\ < 90 deg lines of sight crossing the 
dust plane. Of course, this effect is not observed when one 

2 We derive these luminosities in V with the formula L(M V ) = 
C v Iq-oam v x An(D^ Wpc) f, where C v = 1.97 x lO" 8 W nr 2 . 
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Fig. 3. Plane parallel approximation and cosecant method. The matter (HI gas or stars for instance) is assumed to be organised in planes parallel 
to the Galactic plane, (a): We assume an infinite thick sheet of uniform density centred on the Galactic plane. In this configuration, the column 
density perpendicular to the plane can be derived from the integration of the column density along various directions b = bi, bi or & 3 ... as 
defined in Eq. 5. (b): In this paper, we consider a truncated sheet for the plane parallel configuration, and then study this approximation for an 
exponential disc. 
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Fig. 2. Correction ratio as a function of wavelength necessary to derive 
the surface brightness of the Galaxy (observed face-on) from observa- 
tions at the Solar Neighbourhood. We display the corrections (thick 
lines) for a cylinder of radius 1 kpc and 600 pc (corresponding to the 
1/ sin/? < 6 cut for the cosecant law approximation). We also provide 
for comparison purposes the corresponding curves (thin lines) for a 
cylinder radius of 300 pc. For the wavelength dependence of the ex- 
tinction, we assume Ry = 3.1 and use the prescription of Fitzpatrick 
(1999). The different wavelengths for which we project in cosec/? all- 
sky data are indicated (and numerical values provided in Table E. 1). 



computes the surface brightness for an observer of the face-on 
Galaxy, as only the lines of sight 3 (\b\ = 90deg) perpendicu- 
lar to the Galactic plane are integrated. In Fig. 2, we display 
the correction factor necessary to convert our local estimates 
(measured from the Solar Neighbourhood) into surface bright- 
ness of the face-on Galaxy (see also Tab. E.l). 



2.1.2. Integration in a "cone" 
approximation 



cosecant-law 



We integrate the flux received at the Sun in a cone, according 
to the cosecant-law approximation (see Fig. 3). The slope of 
the corresponding function v^O^I) is the surface brightness of 
our Galaxy defined above, assuming an homogeneous stellar 
density galaxy. 



v v Iv(\b\) 



f n{ 

JrJ.Mv 



2x n(r,bJ,M v ) 

'r,I,Mv 

f(A v (r,b))L(M v ) 



An 



cos bdb dr dl dM\ 



(4) 



3 For our Galaxy, the extinction corresponding to the total disc 
thickness is 0.12 Vmag. 
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Fig. 4. Simulations. We present the cosecant laws computed for expo- 
nential disc and plane parallel geometries. We also consider the effect 
of extinction. When not indicated otherwise, the slopes are measured 
for 1 / sin |fc| < 6. The fixed and variable h- geometries are barely dis- 
tinguishable. 
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Fig. 5. Simulations. V Surface Brightness derived from the cosecant 
laws as a function of the chosen cut in 1 / sin j£>j. The symbols and 
lines follow the same convention as Fig. 1. The line widths and types 
are also defined in Fig. 4. 



The normalisation factor in front of Eq. 4 is due to the nor- 
malisation of this method with respect to the integration in a 
cylinder. As explained in Appendix C, it has been determined 
for a uniform plane parallel geometry. Figure 4 displays the 
vvlv(\b\) function, used to compute the surface brightness. To 
a good approximation, the curves are linear. 

In the case of a plane parallel approximation without ex- 
tinction, the slope of the Iv(b) function with respect to 1 / sin b 
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Table 1. Simulated V surface brightness expressed in different units. 
The first and second columns provide the surface brightness estimates 
derived from the Solar Neighbourhood (S.N.) and for the face-on 
Galaxy observed at the Solar Neighbourhood. The third column pro- 
vides the estimate in the absence of extinction. The last column gives 
the unit of value provided in the different lines. 



from S.N. 


Face-on @ S.N. 


without extinction 


units 


2.91 x irr 6 


4.72 x irr 6 


4.98 x irr 6 


Wm~ 2 


9.7 


15.6 


16.5 


L oP c- 2 


23.90 


23.37 


23.31 


mag arcsec~ 2 



is equal to the value at b = 90 deg, namely: 
I v (b)-I v (b = 90 deg) 



I v (b = 90 deg) 



(5) 



For a linear function, Eq. 5 is an equality and is equivalent to 
vyS v = Iv(b = 90 deg) = Iv(b) sin b, where v v I v (b = 90 deg) 
corresponds to the surface brightness. Figure 4 shows that, 
without extinction, the exponential disc geometry is very close 
to this approximation. In the presence of extinction, the two 
geometries (plane parallel and exponential disc) are also very 
close. Figure 5 displays the sensitivity of the cosecant-law ap- 
proximation to the cut in (1/ sin \b\) M AX- 

2.1 .3. Comparison 

Plane parallel and exponential disc geometries provide very 
close results. In Figure 1, we compare the two modes of in- 
tegration in V. 

In the absence of extinction, the two modes of integration 
(in a cone and in a cylinder) are equivalent for R cy i < 2kpc. 
With extinction, a flattening is observed for both methods. For 
integration in a cone, it is sensitive to the cut in 1 / sin b used 
to compute the slope (see Fig. 5). For integration in a cylinder, 
the surface brightness decreases as the radius of the cylinder 
considered increases (see Fig. 1). 

As displayed in Fig. 1, we estimate that 1 / sin \b\ < 6 is the 
best compromise to sample the surface brightness of the Solar 
Neighbourhood over the whole spectral range. The cosecant- 
law approximation (computed with 1 / sin \b\ < 6) is compati- 
ble with the surface brightness computed in a cylinder of radius 
1 kpc. The plane parallel model provides a similar equivalence 
at 600 pc. For this simple geometry, this radius can be com- 
pared to the maximum radius probed by the integration in a 
cone: h z / tanb mal . For the parameter values considered, stars 
up to a radius of 710pc are probed. For the exponential disc, 
we expect a similar difference (but asymmetric). 

Table 1 provides the values thus obtained for the V surface 
brightness of our Galaxy at the Solar Neighbourhood for the 
different cases. It is important to note that the difference be- 
tween the surface brightness estimated without extinction and 
the face-on Galaxy surface brightness (with extinction) is so 
small that scattered light can be neglected (less than 6% of the 
surface brightness). 
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Fig. 6. Simulations in V with extinction. Integrated V flux in 1 mag 
bins as a function of 1/ sin \b\ for the bright stars (V < 6). 
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Fig. 7. Simulations in V with extinction. Integrated V flux in 1 mag 
bins as a function of 1/ sin \b\ for fainter stars (V > 7). 

2.2. Simulation of a catalogue reduction 
(observational procedure) 

In the following, we will compile several whole sky stel- 
lar catalogues to measure the surface brightness of the Solar 
Neighbourhood. In order to understand the possible trends that 
could affect this observational procedure, we simulate this re- 
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Fig. 8. Simulations in V with extinction. Contribution of each appar- 
ent V magnitude bin to the surface brightness for 4 sky areas. On the 
left (resp. right) panel, the stellar density (Eq. 1) is computed with a 
variable (resp. fixed) scale height h z (M v ) (resp. h z = 120 pc). 



duction for a V catalogue based on the cosecant-law approxi- 
mation. We take our exponential geometry model with extinc- 
tion and a variable scale height hz(My) and decompose the 
I(\b\) function into apparent V magnitude bins. We also divide 
the sky into 4 areas, in order to study asymmetries between the 
North and South Galactic hemispheres and the Galactic Centre 
and Anti-Centre directions. The sum of the slope of each bin 
provides the surface brightness at the Solar Neighbourhood. 

Figures 6 and 7 display this simulation for V < 6 and V > 7 
stars. The flux is higher in the direction of Galactic Centre. We 
observe a significant flattening for large 1 / sin due to ex- 
tinction. We also clearly detect the asymmetric position of the 
Sun above the Galactic plane (z = 15pc). The left (resp. right) 
panel of Fig. 8 summarises the contribution of each apparent 
V magnitude bin to the surface brightness for hz(My) (resp. 
h z = 120 pc). Table A.l summarises the corresponding contri- 
bution per apparent magnitude bins. It also provides the fitted 
value at b = 90°. 

Each bin with V < 2 and V > 21 stars contributes to no 
more than about 1% of the surface brightness, while the main 
contributors (86.5% of the total surface brightness) are 5 < V < 
16 stars. Uncertainties on the very bright stars, which are often 
missed in stellar catalogue (due to their saturation, extinction 
in the stellar disc, etc.), and dim stars (close to the detection 
limit) will not bias the surface brightness estimate. 

Note that this method, uncorrected for extinction, provides 
an estimate of the surface brightness computed from the Solar 
Neighbourhood, and has to be corrected as indicated in Fig. 2 
and Table E. 1 to obtain a face-on value. 
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Fig. 9. Corrected direct and indirect synthetic estimates of the sur- 
face brightness of our Galaxy at the Solar Neighbourhood. In the top 
panel, we display the two estimates obtained with the compilation of 
the V < 7.3 Hipparcos Catalogue stars. For the indirect estimate, we 
indicate the contribution of the different types of stars. In the middle 
panel, we display the contribution of M v > 7 stars with our simula- 
tion tool for different configurations. We also consider ± lcr luminosity 
functions to estimate the scatter due to uncertainties in the luminosity 
function. These dim stars mainly contribute to the near-infrared part of 
the spectrum, but this integrated surface brightness remains an order 
of magnitude below their bright counter-parts. In the bottom panel, 
we study the influence of the extinction, perpendicular to the Galactic 
plane, which might have a marginal influence in the UV. We assume a 
fixed extinction A v applies to the whole population. 



3. Stellar Surface Brightness: a synthetic 
approach with the Hipparcos data 

The Hipparcos Catalogue (Perryman et al. 1997; van Leeuwen 
et al. 1997) with the Hipparcos Input Catalogue (HIC) (Turon 
et al. 1992, 1995) constitute a unique database to explore the 
Solar Neighbourhood. Spectral type, parallax, and V magni- 
tude are provided for each V < 7.3 star. We restrict the sample 
to 13 374 stars with a 3<x parallax and with the constraints on 
the completion limits described in Appendix D. In the follow- 
ing, we use the stellar spectral flux library from Pickles (1998) 
to produce synthetic spectra: a spectrum is associated with each 
star and normalised to the observed V flux. Hipparcos stars are 
then combined in two different ways to get estimates of the 
surface brightness. In Sect. 3.1, we directly sum the flux of 
the observed stars with a sin \b\ weighting. However, this direct 
measurement requires a volume correction due to the apparent 
magnitude limit of the sample used. This corrected direct mea- 
surement corresponds to the cosecant averaging (integration in 
a cone). In Sect. 3.2, we compute the surface brightness with 
the luminosity of the stars with the integration in the cylinder 
approach. Appropriate volume corrections are also required. 
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3.1. Corrected direct measurement of starlight: a 
synthetic spectrum 

The cosecant-law approximation (see Eq. 5) enables to com- 
pute directly the surface brightness as the sum of the stellar 
fluxes weighted by sin \b\. Even though the amount of dust is 
small at the Solar neighbourhood, it significantly affects the 
surface brightness computation (see Fig. 1). We have shown 
in Sect. 2 that the surface brightness of our Galaxy observed 
face-on is very close to the surface brightness computed for a 
Galaxy without dust. Hence, we subtract the extinction to each 
star relying on our dust modelling (see Appendix E), and obtain 
the lower bound: 

all stars 

vS v > 2x vF v (M v ,SpT)/f(A v )xsm\b\ 



all stars 



2 * Z 



L A (M v ,SpT) 



. , x sin \b\ 
Anr 2 



(6) 



where F V (M V , S pT) (resp. L A (M V , S pT)) is the flux (resp. lu- 
minosity) spectrum in W trT 2 Hz -1 (resp. W) corresponding to 
a star with a magnitude My and a spectral type SpT, and r, the 
actual distance of this star. 

According to our simulations (Table A.l), the cut in ap- 
parent magnitude (V < 7.3) limits our sensitivity to 17% of the 
total V stellar surface brightness. We apply a volume correction 
relying on the knowledge of the Hipparcos parallaxes together 
with the well-defined completeness limit (V < 7.3): each abso- 
lute magnitude bin (My - \ < m < My + i) is complete within 
a radius R(My + |) (see Appendix D.2). We compute 



M V <7 X;£ 



L A (xi,SpT)sm\b\ LF 
~ x h 



[My-\,Mv+\[ 



Anr 2 



, w: p \R max ) 
'w: ph (R( Xi )) 



(7) 



here L A (xi, S pT) is the spectrum (in W) corresponding to a 
star with a magnitude My - \ < Xi < My + j and a spec- 
tral type SpT, f t LF is the correction of the luminosity function 
provided in Fig. D.4 applied for stars with My < -1.5, and 
Rmax = lkpc is the cylinder radius discussed in Sect. 2. The last 
factor corresponds to the volume correction with: W" ph (y) = 
f n(r, b, I) sin \b\/r 2 dv, dv = cos(b) db dldr. Wl ph (y) corre- 
sponds to the stellar surface density weighted by sin \b\/r 2 com- 
puted within a sphere of radius y. Wt ph R(xi)) is computed 
for each magnitude bin, while Wl ph (R max ) is obtained with a 
weighted sphere with the radius of the cylinder radius R max . 
The corresponding spectrum is displayed in Figure 9. We check 
that if we take W'l ph (5R max ) instead of W'l ph (R, nax ), the estimate 
of the surface brightness increases by less than 14.5%. 

This measurement takes into account the real distribution 
of stars but assumes an smooth and homogeneous distribution 
of stars of the LF along the exponential profile of the Milky 
Way for the corrections. 



3.2. Indirect measurement of starlight: a synthetic 
spectrum 

The second possibility is to compute the surface brightness as 
the mean stellar luminosity spectrum per unit area at the Solar 
Neighbourhood. We rely on Hipparcos data for intrinsically 
bright stars (My < 7) and perform the appropriate corrections, 
as follows: 



M V <1 



V/v 



I E 

' [M v -±,M v +\[ 



L A (xi,SpT) LF 
nR(xi) 2 Ji 



r cyl 



r sph 



(8) 



where L A (M V ,S pT) is the spectrum (in W) corresponding to 
a star with a magnitude My and a spectral type SpT, f FF is 
the correction of the luminosity function. The last factor corre- 
sponds to the volume correction for each magnitude bin. lf yl is 
the stellar surface density computed for the cylinder of radius 
Rmax, while Y? ph is the stellar surface density integrated over a 
sphere of radius /?(*,-) normalised to the surface nR(Xi) 2 . 

This measurement assumes a smooth homogeneous distri- 
bution of the stars and a uniform LF along the exponential pro- 
file of the Milky Way. 



3.3. Discussion 

Figure 9 presents the corrected direct and indirect estimates. 
The two curves are very close in the optical and near-infrared, 
while they differ significantly in the UV. This difference, which 
cannot be explained with the extinction (see bottom panel of 
Fig. 9), is understood as due to the inhomogeneous distribu- 
tion of the OB associations (de Zeeuw et al. 1999, see also 
Appendix D.3). If we remove the very bright stars (My < 
-4.5), the corrected direct estimate is reduced in the UV. If 
we remove the bright stars (My < -1.5), all the discrepancies 
are removed. 

The corrected direct estimate, based on the fluxes, is very 
sensitive to the anisotropic distribution of OB associations, 
while the indirect estimate is more robust, as it is based on 
luminosities, but does not account for the real spatial distri- 
bution of stars. The uncertainties on the bright end of the LF 
(see Fig. D.l) are also a source of error. By chance, our indi- 
rect estimate is in good agreement with the UV measurements 
of Gondhalekar et al. (1980) (see Fig. 16). However, it is a dif- 
ferent quantity that has been measured: the sum of the UV light 
over the whole sky (S2/68 Sky-survey TD1 -satellite data) with- 
out any extinction correction. This might mean that, due to an 
obvious selection bias, we detect only stars that are moderately 
extinguished, but this might also be a coincidence as the mea- 
sured quantities are different. 

The significant difference of our two estimates in the UV 
range is due to the irregular distribution of OB stars (considered 
homogeneous in the indirect estimate). 

The surface brightness obtained in the optical can be com- 
pared with the study of Flynn et al. (2006). These authors use 
accurate data on the local luminosity function and the disc's 
vertical structure to measure the following surface brightnesses 
(see their Table 5): fiB = 23.51 mag.arcsec -2 , fly = 22.93 
mag.arcsec -2 , and pj = 22.03 mag.arcsec -2 . These values are 
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Table 2. Contribution of the B < 6. stars to the surface brightness (see 
Fig. 14). 2 783 Hipparcos stars are thus been used. Values with the 
observed extinction are provided together with the values corrected 
for extinction. 



A (fim) 



0.440 (Bj) 



With extinction 
5.640 xlO- 7 



Without extinction 
6.186 xlO"' 



in good agreement with our corrected direct estimates (see 
Tables 6 and 7). Our indirect estimates are slightly weaker for 
the reasons discussed above. 



4. Optical surface brightness 

4.1. USNO-A2 and TYCHO-2 data 

The USNO-A2 B and R magnitudes have been calibrated as 
explained in Appendix G. 

4.1 .1 . B measurements 

Figures 10 and 1 1 display the cosecant laws obtained for 4 sky 
areas with the USNO-A2 and TYCHO-2 stellar catalogues and 
for each magnitude bin. Table A. 2 summarises the slopes mea- 
sured with the fit to the cosecant laws and provides a com- 
parison with simulations. These simulations are based on the 
method described in Sect. 2.2. A B luminosity function has 
been derived from the V luminosity function displayed in Fig. 
D.l, assuming that 50% of the V < 2 stars are red giants. We 
relying on Lang (1992, p:137-145) to relate the V absolute 
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-6 -4 -2 
1 /sin (b) 
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-6 -4 
1/sin (b) 



Fig. 10. Cosecant laws for 12 < B < 19 stars from the USNO-A2 
catalogue. Irregularities, significant towards the Galactic centre, are 
minimised by our cut 1 / sin \b\ < 6. (See Table A.2.) 
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Fig. 11. Similar to Fig. 10 for B cosecant laws computed for 6 < B < 
1 1 TYCHO-2 stars. (See Table A.2.) 

Table 3. Contribution of the R < 6. stars to the surface brightness 
computed the same way as those displayed in Fig. 14. 7 010 Hipparcos 
stars are thus considered. Values with the observed extinction are pro- 
vided together with the values corrected for extinction. 




A (irni) With extinction 



0.640 (R) 



6.98 x 10- 7 



Without extinction 
7.41 x 10- 7 



magnitudes to spectral types from Pickles (1998), and to per- 
form this transformation. We thus find excellent agreement be- 
tween our simulations and the TYCHO-2 data. The USNO-A2 
data provide estimates systematically larger than our simula- 
tions. The emission in the Galactic Centre South direction is 
perturbed as displayed in Fig. 10. Note that simulations show 
that B > 19 stars do not contribute more than a few percent to 
the total surface brightness. 

We combine the B USNO-A2 and TYCHO-2 estimates for 
B > 6 stars with the synthetic uncorrected estimate based on 
(B < 6) Hipparcos stars (see Table 2 and Fig. 14), and cor- 
rect the extinction with the global factor provided in Fig. 2 and 
Table E.l, as follows: 

v B S B = (5 .640 x 1 0~ 7 + 3 .004 x 1 0~ 6 ) x 1 .72 
= 6.13 xlO- 6 WtrT 2 = 23.51 magarcsec" 2 



(9) 
(10) 



4.1.2. R measurements 

The data reduction of the R measurements of the calibrated 
USNO-A2 data is presented in Table A.3 and Figure 12. We 
clearly observe an excess of the stellar data with 12 < R < 18 
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Fig. 12. Integrated R fluxes as a function of 1/ sin(fc) computed with Fig. 13. Cosecant laws for Pioneer 10/11 B (full-line) and R (dashed- 
the USNO-A2 data for 12 < R < 18 stars. We decompose the con- line) data (Gordon et al. 1998). This corresponds to V > 6.5 stars and 
tributions originating from different sky areas. No correction has been diffuse emission. There is a hole in the Galactic Anti-Centre North 
applied. (See Table A.3.) area (bottom left) that has been removed. No correction has been ap- 

plied. (See Table 4.) 



with respect to our simulation 4 . This cannot not be accounted 
for the thick disc dominated by dim stars (M v > 4), which 
do not contribute significant to the surface brightness (see Fig. 
D.2). We carefully check with the USNO-B1 (Monet et al. 
2003) catalogue that our calibration is satisfactory. In addition, 
this excess has been confirmed with the CCD UCAC-1 data 
(see Appendix F). 

We complete the 12 < R < 18 USNO-A2 estimate with 
simulations presented in Table A. 3, and correct the extinction 
with the global factor provided in Fig. 2 and Table E. 1, as fol- 
lows: 



Direct measurement computed with the Hipparcos catalogue 




v R S R = (2.205 x 10~ 6 + 2.555 x 10~ 6 ) x 1.41 
= 6.71 x lO^WtrT 2 = 22.6 mag arcsec" 2 



(11) 
(12) 



4.2. Pioneer 10/11 data 

As discussed in Gordon et al. (1998) (see also Leinert et al. 
1998), the Pioneer 10/1 1 data provide an excellent all-sky sur- 
vey in B (0.437/im) and R (0.644yum) for the study of integrated 
flux. In contrast to the Tycho star mapper background analy- 
sis of Wicenec & van Leeuwen (1995), the Pioneer 10/1 1 data 
taken beyond 3.26 AU provide sky-maps devoid of detectable 
zodiacal light. We use here the Pioneer 10/1 1 maps kindly pro- 
vided by K. Gordon (see Gordon et al. 1998 and Witt, Gordon 

4 To account for red giants, we apply in R the modification of the 
luminosity function described in Sect. 4.1.1 to perform this simula- 
tion. 



Fig. 14. Direct contribution of bright stars to the surface brightness. 
These synthetic estimates have been computed with Eq. 6 for V < 6.5 
stars and for B < 6 stars from the Hipparcos catalogue. 8 013 and 
2 783 stars have thus been used, while the few neglected stars do not 
contribute to more than 10% in V (see Table 5). The vertical lines 
indicate the wavelengths discussed in the text (B j and Pioneer B/R). 
Note that the extinction considered here is the extinction affecting each 
individual star not the global extinction perpendicular to the Galactic 
plane. 



& Cohen, in preparation): they correspond to the first iteration 
maps discussed in Gordon et al. (1998). These maps contain all 
relevant flux (starlight and possible diffuse components), ex- 
cept for stars brighter than V = 6.5 (see Toller et al. 1987). 
Figure 13 displays the 1/ sin(fo) laws obtained in B and R for 
the different sky areas, while Table 4 provides the correspond- 
ing slopes. The R fluxes are systematically larger than the B 
fluxes: as discussed in the following, this is a consequence of 
the extinction. 
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Table 4. B(0.437yum) and W(0.644yum) (Pioneer 10/11 derived) surface brightness, given in 10 9 Wm 1 . We give the values fitted for 1/ sin(i) < 
6 (see Fig. 13). The surface brightness thus obtained accounts for the V > 6.5 stars as well as for a possible diffuse component. 



A (jim) 


Galactic Centre 
North South 
\b\ = 90° slope \b\ = 90° slope 


Galactic Anti-Centre 
North South 
\b\ = 90° slope \b\ = 90° slope 


0.437 (B) 
0.644 (R) 


1050. 1148. 
1368. 1684. 


1048. 1444. 
1542. 1958. 

Sum 


1046. 886. 996. 1040. 
1320. 1164. 1500. 1234. 


0.437 (B) 
0.644 (R) 


2098. (\b\ = 90°), 2592. (slope) 
2910. (\b\ = 90°), 3642. (slope) 


2042. (\b\ = 90°), 1926. (slope) 
2820. (\b\ = 90°), 2398. (slope) 




Total B (0.437jum) value: 4140. (\b\ 
Total R (0.644//m) value: 5730. (\b\ 


= 90°), 4518. (slope) 
= 90°), 6040 (slope) 



Table 5. Contribution of the V < 6.5 stars to the surface brightness 
(see Fig. 14). Values with the observed extinction are provided to- 
gether with the values corrected for extinction. We estimate that we 
do not miss 1 more than 3% of the brightness in V. 



A (jim) 


With extinction 


Without extinction 


0.437 (B) 


7.020 xlO- 7 


7.751 xlO" 7 


0.644 (R) 


7.261 xl0~ 7 


7.699 xlO" 7 


V bad stars 1 


2.465xl0~ 8 





1 This corresponds to stars with uncertain parallax or unspecified spec- 
tral type neglected for the previous direct estimates. 



Relying on Eq. 5, we compute the contribution of the bright 
(V < 6.5) stars, independently with a synthetic direct estimate 
based on the Hipparcos Catalogue, which is complete for these 
stars. We nevertheless neglect stars with uncertain parallax or 
undefined spectral type, but estimate that they do not contribute 
more than 10% in V. 

We combine the Pioneer B and R estimates with the (V < 
6.5) stars synthetic uncorrected estimates, and correct the ex- 
tinction with the global factor provided in Fig. 2 and Table E.l 
(see also Appendix E), as follows: 

S B = (7.02 xl(T 7 + 4.52 x 1(T 6 ) x 1.72 

= 8 .98 x 1 0~ 6 W nT 2 = 23 . 1 mag arcsec" 2 (13) 

S R = (7.26 xl0~ 7 + 6.04 x 10~ 6 ) x 1.41 

= 9.54 x lO^WrrT 2 = 22.3 mag arcsec" 2 (14) 

Our estimate is brighter than the estimate of the surface bright- 
ness obtained in B by van der Kruit (1986) with the Pioneer 
data (jiB = 23.8 + O.lmag arcsec~ 2 ), who did not correct for 
extinction. 



The values obtained here are systematic brighter than our 
synthetic estimates and than the local values provided by Flynn 
et al. (2006). As later discussed in Sect. 6, we suggest that this 
disagreement is due to a larger volume probed with the USNO- 
A2 and Pioneer 10/1 1 all-sky surveys. 



x10 Galactic Centre North x 10 Galactic Centre South 




2 4 6 8 10 -10 -8 -6 -4 -2 

1/sin(b) 1/sin(b) 
-4 -4 

x 10 Galactic Anti-Centre North x 10 Galactic Anti-Centre South 




1 /sin (b) 1/sin(b) 

Fig. 15. Cosecant laws for the 1.25 yum COBE/DIRBE data averaged 
over Galactic latitude. These values are derived from Zodi-Subtracted 
Mission Average (ZSMA) data distribution. The sharp peaks corre- 
spond to bright stars (0 < V < 1 and M type). We performed similar 
plots for the other 9 COBE/DIRBE available wavelengths. Wiggles 
observed for 1 / sin |i>| > 6 are due to the COBE/DIRBE sampling 
(pixels of 0.32°x0.32°). 



5. Infrared surface brightness: COBE/DIRBE data 

We consider the Zodi-Subtracted Mission Average (ZSMA) 
maps (Kelsall et al. 1998) from the COBE/DIRBE project. We 
plot the cosecant laws for the wavelengths 1.25, 2.2, 3.5, 4.9, 
12, 25, 60, 100, 140 and 240^m as shown in Figure 15 (for 
1.25/im). We then perform a linear fit (slope and \b\ = 90° val- 
ues) with the parameters provided in Table A.4. We are sensi- 
tive to bright point sources (like Betelgeuse) up to 4.9/im, that 
we eliminate with a 5<x clip when necessary. Unlike Boulanger 
& Perault (1988), we do not subtract nearby dust complexes 
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from the study as we are interested in the light of all the compo- 
nents. Otherwise, we observe quite a linear behaviour except at 
12 and 25 yum where the light emission is very irregular, as ob- 
served by Boulanger & Perault (1988) with IRAS data. These 
data exhibit a large value of the \b\ = 90° intensity and a small 
value of the slope. This large \b\ — 90° intensity is characteris- 
tic of a "non-Galactic" or at least irregular behaviour, and also 
seems to affect marginally the 4.90^111 data. At 12 and 25 yum, 
we observe an unexplained excess for 1 / sin \b\ < 2 (b > 30°), 
which affects the fit. A strange behaviour was also noted by 
Boulanger & Perault (1988) in this area and interpreted as a 
consequence of imperfect Zodiacal light subtraction. Hence, 
we perform the fit for 2 < 1/ sin \b\ < 6 at these wavelengths in 
order to minimise the effect, but these measurements are prob- 
ably still affected by noise and systematic effects. Note that 
Arendt et al. (1998), who studied the Galactic foreground from 
the same ZSMA data, have identified similar spurious effects 
and explained them with residuals of the Zodiacal light sub- 
traction. 

For A < 5/rni, we observe the expected asymmetry between 
the North and South hemispheres due to the position of the 
Sun above the Galactic plane: the integrated radiation inten- 
sity measured in the Galactic Center (resp. Anti-Centre) South 
panel is larger than the one measured in the Galactic Center 
(resp. Anti-Centre) North panel. At longer wavelengths, we ob- 
serve an inverted behaviour: the radiation intensity measured 
in Galactic North panels is larger than the one measured in the 
South panels. This is explained by the presence of a warp of the 
Galactic disc (beyond 12 kpc) towards the North and Galactic 
Centre (see e.g. Burton & Te Lintel Hekkert 1986; Wouterloot, 
Brand, Burton, & Kwee 1990). Last, for A > 60/im, we are 
affected by nearby complexes (that we choose not to remove): 
this tends to enlarge the uncertainties but accounts for local 
fluxes. This probably explains in part the relatively low val- 
ues of the intercepts (\b\ = 90°) found for those wavelengths. 
However, as further discussed in the next section, and displayed 
in Figure 16, the values of the slopes thus obtained are fully 
compatible with those used by Desert et al. (1990) based in part 
on IRAS data, and those computed by Boulanger et al. (1996) 
based on COBE/DIRBE and FIRAS data . 

6. Discussion 

On the one hand, Figure 16 provides a summary of all the 
measurements performed in this paper. We also superpose data 
from the literature as indicated. Perault et al. (1991) provides 
a cosecant measurement in the UV of the diffuse component. 
Interestingly, it confirms that it does not have a substantial con- 
tribution. On the other hand, Table 6 gives numerical values 
of the surface brightness for UBVRI filters and a comparison 
with other works. Table 7 summarises similarly the values for 
near-infrared and infrared. 

We have shown with simulations that the surface brightness 
of the Galaxy at the Solar Neighbourhood as observed from 
outside the Galaxy can be computed with a cosecant averaging 
of all-sky surveys. While face-on estimates performed on exter- 
nal galaxies are not very sensitive to extinction, our measure- 
ments severely suffer from extinction due to our internal view 



point. We have computed the extinction corrections necessary 
to apply to the direct cosecant averaging estimates. Relying on 
these results, we compute synthetic spectral estimates of the 
optical surface brightness. We use the Hipparcos catalogue and 
the Hipparcos Input Catalogue together with a stellar spectral 
flux library to estimate the expected contribution of stars in the 
optical. We use the knowledge of the stellar parallaxes and of 
the completeness limit to perform this estimate. We attempt to 
correct the incompleteness (uncertain parallaxes and stars be- 
low the completion limits) with corrections based on our sim- 
ulations. We then compile optical and infra-red sky-surveys to 
obtain direct observational estimates. 

We observe a systematic difference between our synthetic 
estimates and the optical and infra-red sky-surveys. The es- 
timate based on USNO-A2 R flux in the magnitude range 
12 < R < 18 exhibits a large excess with respect to our sim- 
ulations. This excess is compatible with the disagreement ob- 
served between Pioneer 10/1 1 and our synthetic estimates. On 
the other hand, in the B band the TYCHO and USNO-A2 data 
are in relatively good agreement with the simulations and with 
our synthetic estimate. The Pioneer 10/11 estimate nevertheless 
exhibits an excess (similar to what is observed in R). Various 
arguments could be evoked to provide tentative explanations 
for these differences: 

- The USNO-A2 data exhibit a very irregular pattern in B, 
which is not observed in R. This might be an indication of 
patchy absorption affecting solely the B data. Patchy ex- 
tinction is not included in our simulations nor in the ex- 
tinction correction. However, we would expect the Pioneer 
10/11 data to be affected the same way. There might be 
some other not understood problem with the B data. 

- The synthetic estimates are based on the Hipparcos data 
only (distance limited sample with V < 7.3), while the 
all-sky surveys integrate the whole sky. In practice, the all- 
sky surveys probe an equivalent cylinder radius of 1 kpc, 
while the synthetic estimates are based on the Hipparcos 
data complete within ~1 00-200 pc depending on the mag- 
nitude range. In order to provide a surface brightness esti- 
mate within 1 kpc, we have applied appropriate corrections. 
However, these corrections are normalised to the local stel- 
lar density measured by Hipparcos and assume a uniform 
stellar density following the exponential disc geometry. 
We suggest that this disagreement could be a signature of 
the presence of the local minimum of the stellar density 
linked with Gould's belt. As displayed in Fig. D.3, the 
magnitude interval 10 < V < 18 is roughly composed of 
(My) = 4-5 stars. In addition, we know that the strongest 
contribution comes from -3 < My < 6 (see Fig. D.2). 
The irregular distribution of OB stars (not accounted for in 
our correction for synthetic estimates, nor in our modelling) 
plays an important role. 

- The all-sky surveys estimates follow the same profile as the 
synthetic estimates. This is also compatible with a mini- 
mum of the stellar density (the bright stars are the main 
contributors to the surface brightness). 

- The Pioneer 10/11 data points, when uncorrected for ex- 
tinction, are in good agreement with our synthetic esti- 
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Fig. 16. Summary of the measurements discussed throughout this paper. The data points with no error bars have been computed with the 
cosecant law extrapolation method applied at given wavelengths (COBE/DIRBE, Pioneer 10/11, USNO-A2 & Tycho-2). The Pioneer 10/11 
and COBE/DIRBE 1.25/iin and 2.2yum are corrected for extinction (cf. Table E.l) by the amount indicated by the vertical lines attached to the 
Pioneer 10/1 1 points. In parallel, the dashed and full line curves correspond to the stellar synthetic estimates based on Hipparcos data. In the IR, 
our estimates based on COBE/DIRBE are compatible with Desert et al. (1990), Boulanger et al. (1996) (symbols with error bars). As discussed 
in § 5, the 12/iin and 25/um data points are probably affected by uncertainties due to Zodiacal light residuals. In the UV, the uncertainties are 
large due to the sensitivity to inhomogeneities (OB associations). Our indirect estimate is compatible with the all-sky average performed by 
Gondhalekar et al. (1980). 



mate. However, they seem at odds with respect to the 
DIRBE/COBE J and K data. This is difficult to understand 
as the synthetic estimates are corrected for extinction (star 
by star). 

The recent work of Flynn et al. (2006) based on (partly) 
different data sets and a careful modelling exhibits similar 
trends. Their Tuorla study (which results are quoted in Tab. 

6) probes stars out to circa 200 pc. Their B and V estimates 
are consistent with our corrected direct synthetic estimates, 
since the volume probed with the Hipparcos data is similar. 
In parallel, their I-band estimate is systematically larger than 
our synthetic estimate and closer to our Pioneer 10/11 esti- 
mates. This is in agreement with our Gould's belt interpreta- 
tion as I-band surface brightness should be less sensitive to 
inhomogeneities, and suggests that our incompleteness correc- 
tion, assuming homogeneity of the stellar distribution (see Eq. 

7) , probably misses some stars in I. Our results could explain 
why Flynn et al. (2006) find that the global luminosity of the 
Milky Way appears to be under-luminous with respect to the 



main locus of the Tully-Fisher relation as observed for external 
galaxies. 

We thus conclude that the presence of Gould's belt affects 
our measurements of the surface brightness of the Galaxy at 
the Solar Neighbourdhood. All-sky surveys (Pioneer 10/1 1 and 
DIRBE/COBE) do not suffer from magnitude limits and enable 
to probe the whole volume expected. Our synthetic estimates 
based on Hipparcos data underestimate the surface brightness 
by a factor about 2. We tentatively interpret this difference as a 
signature of Gould's belt. The stellar catalogues tend partially 
to support this interpretation. In the R band, we observed an 
excess in the USNO-A2 data compatible with the Pioneer 10/1 1 
estimate. 
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Table 6. Summary for surface brightness measurements for the usual optical bands and comparison with other estimates. The various 
estimates discussed in this paper have been gathered and provided in the different usual units. As discussed in the text, the values obtained 
for Pioneer 10/11, probing a larger volume than the (corrected estimate) synthetic estimate based on the Hipparcos data, are systematically 
larger (brighter). Our (corrected direct) synthetic estimate is compatible with the Flynn et al. (2006) data in B and V. In the I band, the Flynn 
et al. (2006) estiomate is between our synthetic and Pioneer 10/1 1 estimates. This could be a biais of our incompleteness correction, used 
for synthetic estimates, based on V data, which assumes homogeneity of the stellar distribution. However, the Pioneer data (supported by 
the COBE/DIRBE near-infrared) estimates suggest that the surface brightness of the Milky Way is brighter. This could explain why Flynn et 
al. (2006) observed that the Milky Way appears to be underluminous with respect to the main locus of the Tully-Fisher relation for external 
galaxies. 





Unit 




u 


B 


V 


Rc 


Ry 


Ic 


h 


A 


A 




3660 


4400 


5530 


6400 


6930 


7900 


8785 




WnT 2 




1.48 x 10~ 8 


2.90 x 10~ 8 


1.97 x 10~ 8 


1.44 x 10~ 8 


1.25 x 10~ 8 


9.66 x 10- 9 


7.77 x 1Q- 9 


i4 = 


L 




0.282 


0.558 


0.745 


0.753 


0.737 


0.661 


0.622 










This paper 












Synthetic estimates 




















- corrected direct 


10~ 6 W m 


-2 


3.662 


6.249 


6.081 


5.848 


5.758 


5.46 


5.347 




L Q pc~ 2 




9.08 


15.49 


15.07 


14.50 


14.27 


13.56 


13.25 




Lq PC" 2 




32.20 


27.76 


20.23 


19.26 


19.36 


20.51 


21.30 




mag/" 




23.34 


23.49 


23.10 


22.80 


22.66 


22.44 


22.23 


- indirect 


10~ 6 W m" 


2 


2.433 


4.618 


4.924 


4.864 


4.817 


4.597 


4.495 




L pc~ 2 




6.03 


11.45 


12.21 


12.06 


11.94 


11.40 


11.14 




14 P c ~ 2 




21.38 


20.52 


16.39 


16.36 


16.20 


17.25 


17.91 




mag/" 




23.78 


23.81 


23.33 


23.00 


22.86 


22.63 


22.41 


USNO-A2 


10~ 6 W m" 


2 




6.13 




6.71 










L pc~ 2 






15.19 




16.63 










1-4 pc~ 2 






27.23 




22.08 










mag/" 






23.51 




22.6 








Pioneer 10/11° 


10~ 6 W m 


-2 


(6.23) 


8.98 


(10.34) 


9.54 


(9.79) 


(9.30) 


(9.51) 




L pc~ 2 




(15.44) 


22.26 


(25.63) 


23.65 


(24.27) 


(23.05) 


(23.58) 




14 P c ~ 2 




(54.76) 


39.89 


(34.41) 


(31.41) 


(32.93) 


(34.87) 


(37.90) 




mag/" 




(22.76) 


23.09 


(22.52) 


22.27 


(22.09) 


(21.86) 


(21.60) 


Basic simulations 


10~ 6 W m" 


2 




5.13 


4.23 


4.76 










L pc~ 2 






12.72 


10.51 


11.81 










14 PC~ 2 






22.79 


14.11 


15.68 










mag/" 






23.67 


23.49 


23.02 














Comparison with other works* 










v. der Kruit(1986) 


10~ 6 W m 


-2 




4.68 














L pc~ 2 






11.61 














14 per 2 






20.81 














mag/" 






23.8 












Flynn etal. (2006) 


10~ 6 W m 


-2 




6.12 


7.09 






7.96 






L Q pc~ 2 






15.16 


17.57 






19.74 






14 P c ~ 2 






27.70 


23.59 






29.86 






mag/" 






23.51 


22.93 






22.03 





" Values provided in parenthesis have been obtained as an interpolation based on the corrected direct template adjusted on the Pioneer 
10/1 1 and COBE/DIRBE measurements (factor 1.7). 

* The surface brightness measurements in WirT 2 , L Q pc~ 2 and Lq pc~ 2 provided for the estimates of van der Kruit (1986) and Flynn et al. 
(2006) have been derived with the calibration constants c A and Lq provided on the top of this table and computed with Eq. B.2. These Lg 
pc~ 2 values are compatible with the original values within 10% uncertainties, due to the difference of solar luminosity considered. 
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Table A.l. Simulated V surface brightness (in l(T 9 Wnr 2 ) contributions obtained for each magnitude bin and sky area derived from the 
cosecant laws (see Fig. 6 and 7). We give the values fitted for 1/ sin|£>| < 6. The last column provides the percentages of the total surface 
brightness (PTSB) corresponding to each apparent magnitude bin. 
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Table A.2. Contributions from TYCHO-2 and USNO-A2 stars to the B surface brightness, given in We give the values fitted to 

TYCHO-2 and USNO-A2 B data (see Figs. 10 and 1 1) for each sky area and magnitude bin. We also provide the corresponding values obtained 
with simulations. 
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143.07 


74.60 


175.99 


Total contributions 


Data 


Simulations 




\b\ = 90° 


slope 


\b\ = 


90° 


slope 


-6 < B < 6 








2841.74 


295.43 


6 < B < 11 


2985.77 


978.39 


2577.02 


1106.12 


11 < B < 12 


(98.44) 


(225.75) 


256.04 


251.98 


12 < B < 19 


1167.48 


1476.09 


170.62 


901.59 


6 < B < 19 


4163 




2680. 


5845.42. 


2555.12 


6 < B < 11 and 12 < B < 19 


4153 




2454. 


5845.42. 


2555.12 


6 < B < 19 


(7251 


•) 


(3004.) 


5845.42. 


2555.12 
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Table A.3. USNO-A2 R fluxes, given in r'Wm- 2 . We give the values fitted to USNO-A2 R fluxes (see Fig. 12) for 1/ sin(*) < 6. 
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-0.19 


0.62 
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0.71 


Data sums 


12<fl< 18 


107.79 


701.45 


-84.46 


1108.34 


238.31 


392.02 


254.01 


353.00 


Simulated sums 


-6 < R < 6 


654.35 


99.62 
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638.89 


88.26 


853.29 


103.34 


6 < R < 12 


400.26 


422.13 


499.16 


570.30 


416.02 
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524.58 
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Data 
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\b\ = 
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slope 
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Table A.4. Values of the surface brightness computed with the COBE/DIRBE data, given in KT'WnT 2 . We give the values fitted to the 
COBE/DIRBE data for 1 / sin |*| < 6. The influence of the point sources is removed with a 5<x clipping. The slopes are computed as explained 
for the optical (for 1/ sin \b\ < 6), except for the 12 and 25 yum data. As explained in the text (§ 5), these data exhibit a "non-Galactic behaviour" 
and the corresponding fits are performed for 2 < 1/ sin \b\ < 6. 
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A 
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South 
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Appendix B: Surface brightness definition 



Following the usual conventions, the surface brightness vS v is 
denned as: 



where R is the radius (in m) of the area probed and L A its lu- 
minosity in W at the wavelength A. vS v can be directly trans- 
formed in L^q pc~ 2 . It can also be transformed in a magnitude 
scale jj. (mag.arcsec -2 ) with: 



Ha 



1 X 2 

.51og 10 [vS v x — x ] 



4-n 



(B.2) 



where x, r ^ arcsec 2 = 23.5 x 10~ 12 arcsec 2 /sr is the solid angle 
normalisation. c A is the conversion factor of the observed sur- 
face brightness into magnitude. This conversion factor is for 
instance provided by NED 5 For B, V and R bands, we thus 
have the following conversions: 



Hb = -2.51og 10 [v s S fi x 

//v = -2.51og 10 [vySv x 
Hr = -2.51og 10 [v R S R x 



1 23.5 x 1(T 12 

x : ] (B.3) 



2.90 x 10- 8 
1 

1.97 x 10- 8 
1 

1.44 x 10- 8 



4n 

23.5 x 1(T 12 
4n 

23.5 x 1(T 12 



4n 



] (B.4) 
] (B.5) 



with v b Sb, v v Sy and vrS r the surface brightness in B, V and 
R in WirT 2 , while yU B , ji v and [i R are in B magarcsec -2 , V 
mag arcsec~ 2 and R mag arcsec -2 . 



Appendix C: Normalisation factor of the cosecant 
integration method 

We check in the simple case of a uniform distribution 
(n(r,b,l,M v ) = 1, f(A v (r,b)) = 1 and 4>(M V )L(M V ) = 1) the 
normalisation of the cosecant integration method with respect 
to the integration in a cylinder. 

S',= [ dldzRdR]/nR 2 cvl = 2z max (C.l) 

Jfi=0 Jz=-z„„ Jl=0 



u cosec 



b=-% J 1=0 Jr=i 



2tt r r maI (b) 7 



rdrdl cos b sin \b\db 



4nr 2 

— Zmax (C.2) 

(C3) 

$' x is the simplified version of the surface brightness com- 
puted as the integration in a cylinder, while S' cosec corresponds 
to the integration with the sin \b\ weighting used for the cose- 
cant method. They differ by a factor of 2 understood as the 
difference of the volumes probed. As the behaviour of the ex- 
ponential disc geometry is close to the plane parallel geometry 
(see e.g. Fig. 4), we apply this factor 2 to the normalisation of 
all estimates of the surface brightness based on the cosecant 
approximation, namely the integration in a cone (Eq. 3), the 
corrected direct estimate (Eqs. 6 and 7) and the cosecant com- 
pilations of the various datasets. We check that this factor does 
not change more than 7% with an exponential geometry. 

L A 

V ^ y ~ ~^g2 5 http://nedwww.ipac.caltech.edu/help/photoband.lst 
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Fig. D.l. Luminosity function at the Solar Neighbourhood. The LF 
used for the simulation is presented (Miller & Scalo 1979; Jahreiss & 
Wielen 1997). The LF used for the direct/indirect estimation is also 
displayed with error bars accounting for Poisson noise only. 

Appendix D: Luminosity function (LF) 

We use two estimates of the LF to define <p(My): (1) the Miller 
& Scalo (1979) estimate for My < 2, (2) the Jahreiss & Wielen 
(1997) estimate based on Hipparcos data for My > 2. They are 
displayed on Fig. D.l. We find that this choice of the LF com- 
bined with the variable scale height h z (M v ) of Miller & Scalo 
(1979) is equivalent to a fixed scale height h z = 125 pc (used to 
compare the plane parallel and exponential disc geometries). 

D.1. Uncertainties 

Figure D.2 shows the relative contribution to the surface 
brightness from the Solar Neighbourhood of each absolute- 
magnitude bin. This has been computed with the simulations 
described in Sect. 2.2. Dim stars with My > 7 do not con- 
tribute significantly to the surface brightness (1%), while the 
magnitude interval -6 < My < 6 is the main contributor to the 
surface brightness (98%). It is important to note that the bright 
end (My < 0.5) of the LF constitutes a significant contribution 
(34% of the surface brightness). Hence, theoretical estimations 
of the surface brightness based on the LF depend on the un- 
certainties on the upper tail. However, stars with My < -1.5 
(resp. M v < -3.5) correspond to only 11% (resp. 3%) of the 
surface brightness. This effect disappears when one considers 
stellar catalogues, which provide apparent magnitudes (see Fig. 
8). Figure D.2 displays the value of the absolute magnitude of 
the "mean star" of each apparent magnitude bin, which gives a 
qualitative behaviour of these related parameters. 

D.2. Computation of the LF for direct/indirect estimates 

We present here the estimate of the LF, based on the Hipparcos 
catalogue, used to compute the direct and indirect estimates. 
The Hipparcos Catalogue (Perryman et al. 1997; van Leeuwen 
et al. 1997) is complete for V < 7.3, and parallaxes are avail- 
able for most of the stars. We consider that the number of stars 




Fig. D.2. Contribution of each absolute magnitude bin to the V 
surface brightness (computed with 1 / sin \b\ < 6 from the Solar 
Neighbourhood, with h z (M v )). The upper panel displays the contribu- 
tion to v v S v of each bin of absolute magnitude. The horizontal dashed 
line corresponds to the cumulated value (the estimated surface bright- 
ness). The lower panels exhibit the corresponding decomposition of 
v v lv(\b\) f° r each absolute-magnitude bin (with extinction). These val- 
ues have to be increased by 65% (see Fig. 2) to correspond to a face-on 
value estimate of the V surface brightness. 




Fig. D.3. Mean absolute magnitude stars contribution to each apparent 
magnitude bin. 



in each magnitude bin My - i < m < My + i, where m is 
the absolute magnitude, is complete within a sphere of radius 
R(My + ±) = io(™-( M v+?))/5-+i pc for -6.5 < My < 6.5. We 
count stars in these volumes vol(R) using, 



4>(My) 



N(My -\<m<My + i) 
vol(R) 



(D.l) 



where N(My - \ < m < My + ^) is the number of stars 
found with an absolute magnitude m in this bin, and the vol- 
ume is denned as: vol(R) = J r<R r\(r, b, I; h z (M v ))r 2 d sin bdldr, 
and where rj(r, b, I; h z (My)) = n(r, b, I, M v )/n e (My) (open tri- 
angles in Fig. D.l). This accounts for the exponential profile 
of the Galaxy, with different vertical scales for different mag- 
nitude ranges. We adopt the same scale heights as Miller & 
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A with respect to Miller & Scolo (1979) 
A with respect to Jahreiss & Wielen (1 997) 







Fig. D.4. Fraction f! mssed of the stars missed with respect to Miller 
& Scalo (1979) for each magnitude bin. The correction factor f^ F = 
1/(1 _y™ ss «<) i s applied for the computation of the corrected and indi- 
rect estimates in Sect. 3 in order to be compatible with the simulations. 
The error bars account for the uncertainties on the Jahreiss & Wielen 
(1997); Miller & Scalo (1979) LF used for the simulations and the 
Poisson noise of our estimate. 



Scalo (1979) (h z (M v < 0) = 80pc). We also correct each star 
for interstellar extinction according to our extinction modelling 
described in Appendix E. 

Among the 20906 stars with V < 7.3, we consider the 
18 381 stars with a parallax significant at 3<x. We further re- 
strict this sample by 29% to account the completion limit for 
each magnitude bin. We then treat the 13 374 stars as follows: 

- 9 444 stars have a standard spectral type with a luminosity 
class, that can be found in the Pickles (1998) library. 

- 548 stars either have a standard spectral type with no lumi- 
nosity class available in Pickles (1998) and that has to be 
interpolated between existing templates, or have two spec- 
tral types (the average is then adopted). 

- 3 114 stars have a well defined spectral type, but no lumi- 
nosity class. We rely on My to estimate it, and derive a 
template. 

- 1 Wolf-Rayet star and 20 C/N stars, which do not exceed 
3 .06 x 1 0~ 8 W trr 2 sr" 1 (corrected direct) and 1 .40 x 1 0~ 9 W 
irT 2 sr~' (indirect) in V. 

- 247 stars, with no precise spectral type (63% Am/Ap, 18% 
M), are neglected. We estimate that their contribution to the 
V surface brightness is 6.80 x 10~ 8 W irT 2 sr~' (corrected 
direct) and 5.69 x 10~ 9 W irT 2 sr 1 (indirect). 



Figure D.l summarises the LF thereby obtained. The dis- 
crepancy with respect to Miller & Scalo (1979) is significant 
for luminous stars. The corresponding stellar population is 
strongly inhomogeneous (de Zeeuw et al. 1999), and this part 
of the LF certainly suffers uncertainties larger than Poisson 
noise (see also Appendix D.3). In order, to use the same LF 
as for the simulations, we correct the LF as displayed in Fig. 
D.4, for stars brighter than magnitude -1.5. 

D.3. 3D visualisation of the Hipparcos data 

We present a 3D visualisation of the Hipparcos catalogue on 
http://www.obspm.fr/sbhip06. We also position a schematic 
Gould's belt in the Galactic plane following Guillout et al. 
(1998). Gould's belt is significantly larger than the volume 
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Fig. E.l. A v extinction dependence on Galactic latitude. The extinc- 
tion used for the simulations is compared with other models (Drimmel 
et al. 2003; Hakkila et al. 1997; Mattila 1980b) and the dust maps 
(Schlegel et al. 1998). 



probed by the Hipparcos catalogue. The early-type stars are 
distributed in the Galactic Disk, but do not exhibit a clear 
(20deg) inclination. More generally, the Hipparcos stars distri- 
bution is clearly flatten towards the Galactic Poles, but no clear 
correlation with Gould's belt is detected due to the small vol- 
ume sampled. For the majority of the stars, the cut in magnitude 
is clearly seen as a spheroidal distribution of stars. However, 
the bright stars exhibit a more irregular distribution. The ap- 
parent contours are roughly compatible with the local cavity 
observed in the gas distribution (e.g. Lallement et al. 2003), but 
there is not enough information to conclude about an extinction 
bias. 

Appendix E: Extinction Modelling 

In Figure E.l, we display the dependence of the Ay extinction 
with the Galactic latitude obtained from dust maps by Schlegel 
et al. (1998). We find systematic differences with Hakkila et al. 
(1997), based on a compilation of various models. For the pur- 
pose of this work, we wish to use a simple analytical model to 
describe first order effects linked to extinction. We consider the 
extinction model of Mattila (1980b), based on a 2-components 
model. Compared to the A v extinction derived from Schlegel 
et al. (1998), it tends to overestimate the extinction towards 
the pole (even though this excess remains within the error bars 
quoted by Hakkila et al. 1997), and to underestimate it at inter- 
mediate latitudes. We further tune this type of law, and adopt 
the following empirical model: 

A v (r, b) = a v pl sin \b\ x (-^— ) 3 x (1 - exp(-(r| smb\)lp)) 



b + n 



(E.l) 
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Table E.l. Values of the extinction correction applied to cosecant 
measurements for the main filters used in this paper (see Fig. 2 and 
Table E.l). We consider R cy i = 1 kpc. 



A (jum) 


Correction 


A (jim) 


Correction 


0.44 (B) 


1.72 


0.69 (Ry) 


1.36 


0.55 (V) 


1.52 


1.25 (J) 


1.12 


0.64 (R c ) 


1.41 


2.20 (K) 


1.05 



with a v = 1.48 mag/kpc and /3 = 135 pc. The inverse cube 
power of b has been defined empirically to reproduce the b 
dependence on Fig. E.l. This empirical law is also in good 
agreement with the z-dependence of the extinction modelled 
by Hakkila et al. (1997). 

Numerical values 

The extinction correction used to derive a surface brightness 
measurement of the Galaxy at the Solar Neighbourhood is pro- 
vided in Table E. 1 . 

Appendix F: Comparison of the USNO-A2 with 
UCAC-1 data 
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Fig. El. Galactic plane coverage of the UCAC-1 catalogue (Zacharias 
et al. 2000). We restrict our study to the area indicated with thick lines, 
defined by b < and -90° < / < 0°, and we exclude the area between 
I = -33° and b = -1.1 x 1 + 522°. The intensity corresponds to the 
number of stars with 11 < R v < 16. 



Figure F. 1 displays the area of the sky considered for this 
comparison. Figure F.2 and Table F. 1 compare the 2 data sets 
with simulations. 
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Fig. F.2. Integrated R flux as a function of 1 / s'm(\b\) for the UCAC 
(top) and USNO- A2 (bottom) stars lying in the UCAC- 1 area defined 
in Fig. F.l. These fluxes correspond to the contribution of this area to 
the total flux received at the Sun. No correction has been applied. The 
results of the fits for 1/ sin(|Z?|) < 6 are provided in Table F.l. 

Table F.l. Ru (UCAC-1) and R (USNO-A2) fluxes, given in 
10-"W m 2 sr 1 . We give the values (2) and (3) fitted to R v (UCAC-1) 
and R (USNO-A2) fluxes for 1/ sin \b\ < 6 (see Fig. E.l) over the area 
defined in Fig. F.l, for each magnitude bin (1). Last, we provide the 
ratio (4) between the UCAC-1 and USNO-A2 slopes. No correction is 
applied. 



Magnitude range 


= 90° 


Slope 


\b\ = 90° 


Slope 


(1) 


(2) 


(3) 


(4) 


(5) 




UCAC-1 data 


Ru Simulation 


8 < Ru < 9 


22.49 


8.07 


7.32 


31.20 


9 < Ru < 10 


31.98 


16.38 


5.74 


26.86 


10 < Ru < 11 


26.88 


25.62 


4.22 


22.22 


11 < Ru < 12 


17.56 


31.55 


2.87 


17.50 


12 < Ru < 13 


7.72 


32.98 


1.80 


12.99 


13 < Ru < 14 


0.03 


28.88 


1.03 


9.06 


14 < Ru < 15 


-3.95 


25.61 


0.55 


5.90 


15 < Ru < 16 


-4.87 


16.36 


0.29 


3.59 




USNO-A2 data 


R c Simulation 


11 <R< 12 


13.51 


15.16 


2.87 


17.50 


12 < R < 13 


14.94 


35.80 


1.80 


12.99 


13 < R < 14 


4.00 


39.58 


1.03 


9.06 


14 < R < 15 


-4.45 


34.63 


0.55 


5.90 


15 < R < 16 


-7.39 


29.38 


0.29 


3.59 
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Appendix G: Calibration of the USNO-A2 data 

Given the dispersion of USNO-A2 magnitudes, accounting for 
colour effects would introduce an additional source of noise for 
single wavelength flux measurements. As we are interested in 
integrating the total flux at given wavelengths, we deliberately 
neglect colour effects and perform the calibration between sin- 
gle filters. The red measurements correspond to the plate emul- 
sions 103a-E (North) and IIIa-F (south). The equivalent wave- 
lengths are respectively 6450 A and 6400 A which are very 
close to the Rc filter. The blue measurements (103a-O (POSS- 
I O) for the North, and Illa-J (POSS-I E) for the south) dif- 
fer more significantly with equivalent wavelengths of 4050 A 
and 4680 A. We estimate that a 10% uncertainty is introduced 
when aligning these measurements to the Bj data, which re- 
mains negligible for our purposes. 

We cross-identify the 526280 881 USNO-A2 stars (Monet 
et al. 1998) with those of the GSPC2.1 catalogue (Bucciarelli 
et al. 2001) composed of 305 017 stars with CCD (B)VR pho- 
tometry within a 5" x 5" window. In the USNO-A2 catalog, 
the sky is partitioned into 24 zones of South Polar Distance 
(SPD), each of width 7°.5. We thereby find 257 352 stars spread 
over 21 zones. For each zone, we plot the residuals in RA and 
DEC, and keep the 186568 stars within l<x of the mean of 
each zone. We then compute the relationships between different 
magnitude systems as shown in Figure G.l. The correspond- 
ing coefficients are provided in Table G.l. Busno and /?usno 
magnitudes are thus converted into B j (A„ =4400 A) and R c 
(Ag =6400 A) magnitudes, which are in turn converted into 
fluxes with F BM = C b ,r10- QAxB ' r (C b = 4.26 x 10~ 23 W rrT 2 
Hz~\ C R = 3.08 x 10~ 23 W irT 2 Hz" 1 ). 

Note that we do not detect any significant bias between 
northern and southern hemispheres. Systematic effects are ob- 
served for Rusno data for /?usno > 17.5 and 6 < -14°. They 
will tend to overestimate the contribution of faint stars and are 
most probably an effect of crowding. 
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Table G.l. Coefficients used for the calibration of the USNO magnitudes. In blue (B), B, = a + B(B vmo - 15) + y(5usiM0 - 15) 2 . In red (R), 
Rc = a + B(R\]sno - 15) + yCRusNO _ 15) 2 . The columns are defined as follows: (1) the zone number, (2) the declination band covered, (3) the 
filter used, (4) a, (5) B, (6) y, (7) the RMS error of the polynomial fit, (8) the number of stars effectively used in the fit. We thus used 104 354 
stars in blue and 154491 stars in red. 



Zone 


DEC 


Filters 


a 


B 

H 


y 

1 


Residuals 


N stars 


075 


-82.5° to -75.0° 


B 


15.24 


1.14 


0.0 


0.38 


1111 






R 


15.06 


1.02 


6.99 x 10~ 2 


0.34 


724 


150 


-75.0° to -67.5° 


B 


15.59 


0.95 


2.32 x 10~ 2 


0.45 


2048 






R 


15.05 


1.01 


7.68 x 10~' 


0.42 


1320 


225 


-67.5° to -60.0° 


B 


15.80 


1.00 


-1.67 x 10~' 


0.67 


2006 






R 


14.92 


0.97 


6.45 x 10~ 2 


0.41 


1653 


300 


-60.0° to -52.5° 


B 


15.22 


1.16 


0.0 


0.42 


2509 






R 


15.04 


0.95 


8.09 x 10~ 2 


0.36 


2204 


375 


-52.5° to -45.0° 


B 


15.48 


1.09 


0.0 


0.40 


3488 






R 


15.03 


0.99 


7.31 x 10~ 2 


0.38 


1949 


450 


-45.0° to -37.5° 


B 


15.27 


1.15 


0.0 


0.36 


2206 






R 


15.05 


0.98 


5.13 x 10~ 2 


0.33 


1220 


525 


-37.5° to -30.0° 


B 


15.47 


1.03 


9.07 x 10~ 3 


0.50 


7590 






R 


14.90 


0.88 


8.26 x 10~ 3 


0.47 


5061 


600 


-30.0° to -22.5° 


B 


15.63 


1.00 


9.59 x 10~ 3 


0.48 


8879 






R 


15.01 


0.94 


4.88 x 10~ 2 


0.40 


5771 


675 


-22.5° to -15.0° 


B 


15.58 


1.00 


1.44 x 10~ 2 


0.42 


6065 






R 


14.99 


0.95 


3.96 x 10~ 2 


0.42 


4578 


750 


-15.0° to -7.5° 


B 


15.60 


1.05 


0.0 


0.38 


4196 






R 


15.16 


0.92 


1.21 x 10~ 2 


0.41 


6624 


825 


-7.5° to 0.0° 


B 


15.39 


1.06 


0.0 


0.33 


2470 






R 


15.07 


0.97 


0.0 


0.34 


3949 


900 


0.0° to 7.5° 


B 


15.39 


1.07 


0.0 


0.30 


2171 






R 


15.09 


0.99 


0.0 


0.33 


5047 


975 


7.5° to 15.0° 


B 


15.68 


1.00 


0.0 


0.34 


8428 






R 


15.13 


0.97 


0.0 


0.33 


14084 


1050 


15.0° to 22.5° 


B 


15.47 


0.94 


2.69 x 10~ 2 


0.38 


5945 






R 


15.11 


0.95 


1.68 x 10~ 2 


0.39 


11035 


1125 


22.5° to 30.0° 


B 


15.28 


1.05 


9.92 x 10~ 3 


0.33 


7549 






R 


15.11 


0.96 


1.47 x 10~ 2 


0.33 


12974 


1200 


30.0° to 37.5° 


B 


15.21 


1.13 


-6.01 x 10~ 3 


0.28 


12011 






R 


15.11 


0.96 


8.08 x 10~ 3 


0.33 


17769 


1275 


37 5° to 45 0° 


B 


15.40 


0.99 


1 68 x 1 0~ 2 


0.33 


12665 






R 


15.02 


0.94 


1.37 X 10~ 2 


0.39 


21022 


1350 


45.0° to 52.5° 


B 


15.34 


1.08 


0.0 


0.32 


9248 






R 


15.07 


0.97 


4.56 x 10- 3 


0.39 


17967 


1425 


52.5° to 60.0° 


B 


15.65 


1.08 


-3.69 x 10~ 2 


0.76 


3720 






R 


15.12 


0.93 


0.0 


0.55 


8983 


1500 


60.0° to 67.5° 


B 


15.23 


1.15 


-1.08 x 10- 2 


0.31 


6468 






R 


15.10 


0.95 


0.0 


0.43 


9996 


1575 


67.5° to 75.0° 


B 


15.65 


1.16 


-2.03 x 10- 2 


0.25 


412 






R 


14.61 


1.00 


2.90 x 10~ 2 


0.56 


1031 
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Fig. G.l. Each panel presents the Bj versus B US no and R c versus Rusno relations for each South Polar Distance (SPD) interval corresponding 
to a zone. The polynomial adjustments are superimposed for each relation, and the corresponding coefficients are provided in Table E.l. 
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Fig. G.l. Continued 



